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We study numerically and analytically the properties of a one-dimensional array of parallel Joseph- 
son junctions in which every alternate junction is a vr junction. In the ground state of the array, each 
cell contains spontaneous magnetic flux $ < $o/2 which shows antiferromagnetic ordering along the 
array. We find that an externally introduced 27r-fluxon "I>o in such an array is unstable and fraction- 
alizes into two tt fluxons of magnitude |$o- We attribute this fractionalization to the degeneracy of 
the ground state of the array. The magnitude of the flux in the fractional fluxons can be controlled 
by changing the critical current of the n junctions relative to the junctions. In the presence of an 
external current, the fluxon lattice in the antiferromagnetic ground state can be depinned. We also 
observe a novel resonant structure in the V-I characteristics above the depinning current due to the 
interaction between the fluxon lattice and the array. 

PACS numbers: 74.81. Fa, 75.10.Pq 



One of the exciting developments in the field of 
Josephson devices is the fabrication of the three terminal 
controllable Josephson junction The supercurrent 
through a Josephson junction is given by / = /c sin( A0) 
where A4> is the gauge-invariant phase difference be- 
tween the superconductors, and the critical current Ic 
depends upon the junction geometry, normal state re- 
sistance Rn and the temperature T. Morpurgo et al.^ 
showed that the supercurrent through a superconductor- 
metal-superconductor (SNS) junction changes on pass- 
ing a control current through the normal metal. For 
such a junction, the supercurrent / oc shi(A(f) + x) where 
the additional phase difference % is dependent on the 
current through the normal metal. Further theoretical 
work showed that in the diffusive limit of the junction, 
the additional phase factor x can be made tt, thus re- 
versing the direction of the supercurrent with respect to 
the phase difference Ai^ Josephson junction with 

X = TT is referred as the tt junction (we use the term 
junction for the Josephson junction for which x = 0)- 
The TT junction has now been realized in several experi- 
ments The fabrication of such tunable junctions 
has opened immense possibilities for new applications, 
as demonstrated recently by the development of control- 
lable TT-SQUID 121 . 

The next natural step in this field would be to consider 
Josephson junction array (JJA) containing tt junctions. 
Theoretically, Kusmartsev considered a loop contain- 
ing an odd number of tt junctions and showed that the 
loop contains spontaneous magnetic flux in the ground 
state. In the continuum limit, the long Josephson junc- 
tion with alternating critical current density have been 
studied which shows self-generated magnetic fli ix P , . 
Recent studies of JJAs with tt junctions 0, Il2l| have 
shown some novel features arising out of the interplay 
between and tt junctions. Moreover, JJA is a unique 
system which provides experimental realizations of sev- 



eral interesting physical phenomena, some examples of 
which are field induced superconductor to insulator tran- 
sition Aharonov-Casher effect coherent emis- 
sion of radiation . One is then led to ask as what new 
physical phenomenon exist in the ID JJA containing tt 
junctions. 

In this paper, we study numerically and analytically a 
new class of ID JJA: an array of parallel Josephson junc- 
tions in which every alternate junction is a tt junction. 
The ground state contains spontaneous magnetic flux in 
each cell and are ordered antiferromagnetically along the 
array. We find that a quantum of flux (fluxon) with a 
2tt kink in the phase is unstable in such an array, and 
fractionalizes into two spatially separated tt kink fluxons. 
We also calculate the V-I characteristics which shows a 
novel structure due to resonant interaction between the 
moving antiferromagnetic fluxon lattice and the linear 
waves emitted by the array. 

Consider a ID array of parallel Josephson junctions 
containing alternate tt- and junctions (see inset Fig. 
1). The Hamiltonian for this system is 
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where is the gauge-invariant phase difference across 
the i-th junction. The periodic boundary condition is 
imposed at the two ends of the array such that 0o = 4'n 
{N is assumed to be even). In Eq.(l), the first term 
represents the charging energy and the second term is the 
energy of the induced magnetic field due to finite self- 
inductance of the cell (the effect of mutual inductance 
between the cells is neglected). The last term represents 
the energy associated with the Josephson currents. The 
prefactor for the cos (jn term alternates in sign for odd 
(tt-) and even (0-) junctions. The Josephson couphng 
energy Ej = Ic^q/2tt where Ic is the critical current of 
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The time t is in the units of inverse 

-1 



a single junction. 

plasma frequency = y'^j^, where C is the averaged 
capacitance per unit area of the junction. The effective 
Josephson penetration depth is given by Xj = ( 9^^^°/ 
where Lq is the self inductance of a single cell. The Xj 
determines the screening strength of the array, and is 
related to the SQUID parameter P]^ = Xj^. 
From Eq.(l), the equation of motion for (pi is 



l)'sin( 



-7 = A^,( 



(2) 

where a dissipative term is also added |l6l] . The 

coefficient a — Pc^^^, where /3c — 27ri?^/cC/$o is the 
McCumber parameter. The parameter ^ = lext/Ic rep- 
resents the external current through the junction. For 
the numerical simulation, Eq.(2) is integrated using the 
fourth order predictor-corrector method. The consis- 
tency of the steady state solutions was checked using 
different initial configurations of 4>i 's. The magnetic flux 
in the z"* cell is defined as 27r^ = —(0^+1 — <f)i). We 
remark that for the case where all junctions are junc- 
tions (henceforth referred to as the 0-JJA), Eq.(2) is the 
discrete perturbed sine-Gordon equation, and have been 
studied extensively [tR H^. First, we consider the re- 
sults from the numerical simulation. 



we call this array the antiferromagnetic JJA (AFJJA). 
The AF ordering of (and hence, (pi) implies that the 
self-induced screening currents in neighboring cells are 
oppositely oriented. The magnitude of the fiux $ in a 
cell depends on the screening strength Xj as shown in 
the Fig. 1 (inset). With increasing Aj, the magnetic flux 
in the neighboring cells tends to overlap, and $ ^ as 
Xj 00. In the strong screening limit, Xj and 

Next, we consider the consequence of introducing an 
external fluxon in afjja[i3. In the 0-JJA, a fluxon 
corresponds to a 27r-kink in the phase profile ^(a;) and 
the magnetic field (oc ^) is spatially localized on the 
length scale Xj. Figure 2(a) shows the steady state pro- 
files of 27r$i / $0 and (pi in AFJJA in the presence of a 27r- 
fluxon. The self-induced magnetic field of the AF ground 
state has been subtracted from 27r$i/$o- We find that 
a 27r-fluxon is unstable in the AFJJA and fractionalizes 
into two spatially separated fluxons, each carrying half 
the quantum of flux. Also, each of the fractional fluxons 
is a TT-kink in (pi. The magnetic field around the frac- 
tional fluxon decays as exp(— a;/Ae//) where Ae// ~2Aj. 
This should be compared with the 27r-fluxon in the 0- 
JJA [Fig. 2(b)] where Ae// « Xj. The increase in Ae// 
in AFJJA is a consequence of the magnetic flux in the 
AF ground state. 
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FIG. 1: The self-induced magnetic flux 27r$i/<I>o along the 
array in the ground state. The Aj — 1.0 and A'' — 100 (only 
half the array is shown for clarity) . The right inset shows the 
array geometry and the arrows represents the antiferromag- 
netic ordering of the magnetic flux indunced in the cell. The 
left inset shows the dependence of 27r-'^ on A,/. 



Figure 1 shows the ground state flux conflguration 
27r$,/$o for N = 100 and Aj = 1.0. The self-induced 
magnetic flux ^i changes sign across neighboring cells 
with \^i\ — $ remaining constant. Such a conflguration 
of $i is reminiscent of the ground state in ID classical 
Ising model with antiferromagnetic (AF) coupling, as de- 
picted schematically in the inset of Fig. 1. Therefore, 
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FIG. 2; (a) The magnetic flux and the average phase 

profile along the array in (a) AFJJA and (b) 0-JJA, in the 
presence of a 27r-fluxon. The background AF ground state 
has been subtracted from 27r5^ in (a). The schematic in 
inset (a) shows how the fractional fluxons (thick arrows) in- 
terpolate between the two degenerate ground state of $i's. 



It is possible to vary the magnitude of the magnetic 
flux in each fraction by changing i* — 1^ /Ic where 
and are the critical currents of the tt and junctions. 



3 



respectively 12^. Figure 3(a) shows the spatial profile of 
the fractional fluxons for i* — 0.8 and Aj = 1.0. The 
phase change across the fractional fluxons is not tt but 
is dependent on the value of i*. The total phase change 
across both the fractions is always 27r, as required by the 
flux conservation. Figure 3(b) shows the magnitude of 
the integrated flux 271^ in each fractional fluxons as a 
function of i* for Aj = 1.0. The fractionalization occurs 
for i*i < i* < where i*i and i*^ are the two criti- 
cal values. The slope at the critical values «*j and 
j*„ appears to diverge, suggesting a transition between 
the fractionalized state and the single fluxon state. In 
experiments, the magnitude of the flux at the center of 
the fluxon can be measured more easily. Figure 3(b) 
shows the behavior of $m(i*). 

In Fig. 4, we show the numerically obtained param- 
eter space Xj-i*. We have assumed that i* can be var- 
ied independent of Aj. The region of fractional fluxons 
is bounded by i^iiXj) and j*,u(Aj). It is easy to under- 
stand the absence of fractional fluxons in the limit i*—>-0 
since the array becomes a 0-JJA (with lattice constant 
twice the original array) which allows only 27r-fluxons. 
In the opposite limit 1*^00 such that 7° ^ and 
is finite, there are two tt junctions in each cell and the 
array can be shown to be equivalent to the 0-JJA, and 
the fractionalization is again not expected. In obtaining 
the parameter space in Fig. 4, is assumed to be finite 
and fixed which leads to fractionalization for Aj < 0.7 
even as i* 00. 

The simulation results discussed above can also be un- 
derstood analytically. Consider the case of i* = 1 . Define 



l)2m+l 



and 



02m 



(3) 



where m = 0, ^ ~ 1- Thus, Um and Vm are the gauge 
invariant phase differences across the tt and junctions, 
respectively. In the absence of any external fluxon, the 
Um and Vm are invariant on translation by the lattice 
vector along the array. Hence, substituting Um = u and 
Vm — V, Eq.(2) becomes 



sint; = 2Aj(u — v). 



(4) 



There are two trivial solutions of the Eq. (4) : u = v = Q 
and u — v^TT. The non-trivial solution of Eq.(4) is given 
by 



and 



2\^j{tt - 2v). (5) 



For a given value of Aj, the quantities v and u can be 
calculated graphically from Eq.(5). It can be easily ver- 
ifled that the non-trivial solution is the ground state for 
any finite Aj. The magnetic flux in the cell is given by 
$ = ±(u — v)^, where the -I- and — sign is for the 
cell to the left and the right of the tt junction, respec- 
tively. Thus, the magnetic flux alternates in sign along 
the array. The values of u, v and j^j obtained from 



0.5 



(a) 100, 1. = 1.0 

i*= 0.8 



<5 



Q I I I 1^- \ L. 

40 / 




FIG. 3: (a) The spatial profile of the fractionalized external 
fluxon for i* = 0.8 and Aj = 1.0. The magnitude of flux in 
the two fractions are not equal, (b) The integrated total flux 
2n-^ in the two fractional fluxons (represented by dashed 
and dotted lines) as a function of i* = I^/I^. Also shown is 
the maximum flux 27r5r^ at the center of the two fractional 
fluxons $a and $b (represented by symbols). 



Eq.(5) are in excellent agreement with the numerically 
obtained values. In the strong screening limit Aj — > 0, 
u — TT and V — 0, and the flux in each cell attains the 
maximum value ^o/2. In the limit Aj — > 00, the solu- 
tion of Eq.(5) is u — v= Tr/2 which is degenerate to the 
trivial solutions of Eq. (4) , and the ground state contains 
no spontaneous magnetic flux. 

To understand the fractionalization of a 27r-fluxon, we 
note that the AF ground state solution of Eq. (4) is two- 
fold degenerate. If (/> = {u,v} obtained from Eq.(5) is 
one solution, the other solution is obtained by translat- 
ing (f) by one lattice constant. Thus, the other solution 
is (j)' — where u' — v + n and v' = u + t:. The 

degeneracy of the ground state has an important impli- 
cation: the elementary excitation for the array is a kink 
(or domain wall) in which interpolates between the 
two degenerate ground states (f) and 0'. It can be eas- 
ily verifled that the phase change across the 'kink' is tt 
and corresponds to an additional flux $0/2 in the array. 
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FIG. 4: The parameter space Xj-i*- The circles are the nu- 
merically obtained values whereas the full lines are the ana- 
lytical result. 



In the presence of a 2t: fluxon, the energy minimization 
leads to two 7r-kinks which are separated by the degen- 
erate ground states. This is shown schematically in the 
inset of Fig. 2(a). Fractionahzation of a 2tt fluxon in 
AFJJA is thus a consequence of the degeneracy of the 
ground state. A similar phenomenon is observed in poly- 
acetylene (2lJ and certain field theories 22, 2^ where 
fractional topological excitations occur and are related 
to the ground state degeneracy. 

Equations (4) and (5) can be extended for the case 
i* ^ 1. We find that the AF ground state is stable 
and fluxon fractionahzation occurs for i* > 1 when 
Aj < Xju = ^(1 - and for < 1 when 

Aj<Aji = ^(Y^)i/2. This is plotted in Fig. 4 and is in 
good agreement with the numerically obtained behavior 
of iciXAj) and z*;(Aj). Further details of the analytical 
calculations will be given elsewhere. 

Next, we study the dynamical properties of the 
AFJJA. We restrict the analysis to the case z* = 1.0 and 
in the absence of any external fluxon. The V-I curve is 
obtained by sweeping 7 in small steps, and calculating 
V — a{^) in the steady state {V is in units of Rnic)- Re- 
call that for the 0-JJA, all junctions switch from the su- 
perconducting state to the normal state simultaneously 
at 7 = 1. In the AFJJA, the magnetic flux in the ground 
state alters this behavior significantly. 

Figure 5 shows the V-I curve for iV = 30 and Aj = 1. 
The a = 0.1 for the rest of the discussion below(2Jj. The 
important feature of the V-I curve is the appearance of 
a plateau in V above a depinning current 7c. The tran- 
sition to the normal state occurs at a higher current 7,1 . 
By analyzing the spatio-temporal dynamics of fluxons 
on the voltage plateau, we find that the interpenetrat- 
ing lattice of the flux ($) and the antiflux (— <&) move 
in the opposite directions which appears as a stationary 
wave of breathing flux-antiflux pairs |^ . This is evident 
from the 2D space-time plot of 27r^ which is shown for 



a section of the array in the inset of Fig. 5. The ^ in 
each cell oscillates between the positive and the negative 
values, and is in antiphase with the neighboring cells. 
Thus, at any instant of time, the total flux in the array 
is zero as expected from the ground state. We also find 
that for Aj > 1, a linear flux flow regime appears before 
the plateau, whereas for small Aj, the V shows a sharp 
step as shown in the inset of Fig. 5. The Aj dependence 
of the width of the voltage plateau {A'j)s = In — Ic is 
shown in Fig. 6(a). The (A7)s is non-monotonic and is 
maximum for Aj w 2. Below a cut-off A} ~ 0.45, the 
plateau in V disappears and ^c — In- 
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FIG. 5: The full V-I curve of the AFJJA. The direction of 
the current ramp is indicated by the arrow. Upper inset: 
the V-I curves for different values of Aj (the V branch with 
increasing 7 is shown). Left inset: the 2D space-time plot of 
27r-2^ in six cells of the array for 7 = 0.3, where the maximum 
value is 2.2 (black) and the minimum value is -2.2 (white). 

The origin of the plateau (or step) in V can be under- 
stood from the fluxon dynamics. In a 0-JJA, the mo- 
tion of a single 27r-fluxon leads to the emission of small 
amplitude linear waves (plasma waves) due to the dis- 
creteness of the array. The resonances between these lin- 
ear waves and the periodic motion of the fluxon causes 
a series of plateaus in V 0|. Extending this to the 
case of AFJJA, the plateau in V-I can be attributed to 
the phase locking between the moving fluxon lattice and 
the linear waves emitted by the array. The frequency 
ujs of the linear waves can be calculated from Eq. (2). 
In the absence of any external fluxon, the symmetry of 
the ground state allows only waves with the wavevec- 
tor k — 2-K/(2a) to be coupled resonantly to the moving 
fluxon lattice (here a is the lattice constant). Thus, all 
and TT junctions oscillate with the same amplitude and 
phase. Linearizing Eq.(2) for the tt and the junctions 
using u ~ uq eyi^iiuj st) and v — VQe:iq){iujst), respec- 
tively. 



- oj'^u — u = 2Aj(w — u), 
-uj'^v + V = 2Aj(u — w). 



(6) 
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For simplicity, we have used a = 7 = 0. Adding the 
above equations gives the relation between u and v, u = 
v{l — uj'^)/{l + Wg). Substituting for u in the second 
equation leads to 

{l-uj',){l+ujf}+4X'jUj',=0. (7) 

The above quadratic equation in oj'^ can be solved to 
obtain the frequency of the linear waves 




The Us is in units of the plasma frequency ujp. The 
condition for phase locking of the linear waves with the 
moving fluxon lattice then becomes uigT — 2tt, where T 
is the time period corresponding to the motion of the 
fluxon lattice. From the simulation, the T (and hence 
LUs) on the plateau can be obtained from the time evo- 
lution of the magnetic flux in a cell. In Fig. 6(b), 
the LOs calculated analytically is compared with that ob- 
tained from the simulation for increasing Xj. A rea- 
sonable agreement is observed over a range of Aj. The 
discrepancy appears as Xj approaches X*j below which 
the plateau in V is not observed. The voltage Vg on the 
plateau is given by Vs — aus- This is shown in Fig. 6(c) 
and is in good agreement with the simulation. 





FIG. 6: (a) The Xj dependence of the width of the voltage 
plateau (A7)s. (b) Resonance frequency ujs{\j), and (c) the 
voltage Va{Xj). (d) The depinning current 7c. The open 
symbols are from the simulation whereas the bold lines in 
(b), (c) and (d) are the analytical results described in the 
text. 

Finally, the dependence of 7c on the screening strength 
is shown in Fig. 6(d). A finite 7c in JJA is attributed to 
the pinning potential created at the center of the cell . 
For AFJJA, the pinning potential or the energy barrier 
can be defined as AE — Em — Ea, where Em is the 
energy of the array with the fluxon lattice placed on the 
junctions and Ea is the ground state energy. Equating 



the pinning force 2AE to the Lorentz force required to 
overcome the energy barrier gives the depinning current 
7c (per junction), 

j, = l(^2cosiv)-^X^j{7:-2v)^y (9) 

The phase v in the above equation is the solution of the 
Eq.(5). Fig. 6(d) compares the above expression with 
the 7c obtained from the simulation. The agreement is 
good at large Aj, and the deviation appears as Aj A} 
below which no flux flow is observed. 

In conclusion, we have introduced a new class of 
JJA containing tt junctions. We considered the one- 
dimensional case in which every alternate junction is a tt 
junction, and studied its properties numerically and an- 
alytically. The ground state of the array contains spon- 
taneous magnetic flux in each cell and are ordered anti- 
ferromagnetically along the array. A 27r-fluxon in such 
an array is unstable and fractionalizes into two spatially 
separated 7r-fluxons. The fractionalization is related to 
the ground state degeneracy. The V-I curve shows a 
voltage plateau due to resonant interaction between the 
linear waves emitted by the array and the moving anti- 
ferromagnetic fluxon lattice present in the ground state. 
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Note : After the completion of the work, we became 
aware of the experimental fabrication of an array of 
coupled TT-loops using YBaCuO-Nb zigzag structure in 
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